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Abstract 

A quantum compiling algorithm is an algorithm for decomposing ("compiling") an 
arbitrary unitary matrix into a sequence of elementary operations (SEO). Suppose 
Uin is an A^^^-bit unstructured unitary matrix (a unitary matrix with no special sym- 
metries) that we wish to compile. For Nb > 10, expressing SEO requires 
more than a million CNOTs. This calls for a method for finding a unitary matrix 
that: (l)approximates Uin well, and (2 is expressible with fewer CNOTs than Uin- 
The purpose of this paper is to propose one such approximation method. Various 
quantum compiling algorithms have been proposed in the literature that decompose 
Uin into a sequence of ?7(2)-multiplexors, each of which is then decomposed into a 
SEO. Our strategy for approximating Uin is to approximate these intermediate U{2)- 
multiplexors. In this paper, we will show how one can approximate a f/(2)-multiplexor 
by another t/(2)-multiplexor that is expressible with fewer CNOTs. 
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1 Introduction 



In quantum computing, elementary operations are operations that act on only a few 
(usually one or two) qubits. For example, CNOTs and one-qubit rotations are ele- 
mentary operations. A quantum compiling algorithm is an algorithm for decomposing 
("compiling") an arbitrary unitary matrix into a sequence of elementary operations 
(SEO). A quantum compiler is a software program that implements a quantum com- 
piling algorithm. 

One measure of the inefficiency of a quantum compiler is the number of CNOTs 
it uses to express an unstructured unitary matrix (a unitary matrix with no special 
symmetries). We will henceforth refer to this number as Ncnot- Although good 
quantum compilers will also require optimizations that deal with structured matrices, 
unstructured matrices are certainly an important case worthy of attention. Minimiz- 
ing the number of CNOTs is a reasonable goal, since a CNOT operation (or any 
2-qubit interaction used as a CNOT surrogate) is expected to take more time to per- 
form and to introduce more environmental noise into the quantum computer than a 
one-qubit rotation. Ref.[r proved that for matrices of dimension 2^^ [Nb = number 
of bits), Ncnot > :^(4^^ — ^Nb — 1). This lower bound is achieved for A''^ = 2 by 
the 3 CNOT circuits first proposed in Ref.[2]. It is not known whether this bound 
can be achieved for Nb > 3. 
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Suppose Uin is an A^-bit unstructured unitary matrix that we wish to compile. 
As the above table illustrates, compiling Uin is hopeless for Nb > 10 unless we approx- 
imate Uin- We need a method for finding a unitary matrix that: (1) approximates Un 
well, and (2) is expressible with fewer CNOTs than Un- The purpose of this paper is 
to propose one such approximation method. The use of approximations in quantum 
compiling dates back to the earliest papers in the field. For example, Refs.p"^ and 
|3] contain discussions on this issue. As Ref.jS] points out, even when compiling a 
highly structured matrix like the Discrete Fourier Transform matrix, some gates that 
contribute negligibly to its exact SEO representation can be omitted with impunity. 

Refs. [5] [6^ and (7| discuss a quantum compiling algorithm that decomposes 
an arbitrary unitary matrix into a sequence of U (2)-multiplexors, each of which is 



2 



then decomposed into a SEO. Other workers have proposed [H] alternative quantum 
compihng algorithms that also generate ?7(2)-multiplexors as an intermediate step. 

The strategy proposed in this paper for approximating Um is to approximate 
the intermediate f/(2)-multiplexors whose product equals f/j„. In this paper, we 
will show how one can approximate a U (2)-multiplexor by another [/( 2) -multiplexor 
(the "approximant" ) that has fewer controls, and, therefore, is expressible with fewer 
CNOTs. We will call the reduction in the number of control bits the bit deficit 
Sb- Figgis emblematic of our approach. It shows a t/(2)-multiplexor with 3 controls 
being approximated by either a U (2)-multiplexor with 2 controls or one with 1 control. 




Figure 1: Approximating a t/(2)-multiplexor by another [/(2)-multiplexor with Sb 
fewer controls. 



2 Notation 

In this section, we will define some notation that is used throughout this paper. For 
additional information about our notation, we recommend that the reader consult 
Ref. jlUj. Ref.jTU] is a review article, written by the author of this paper, which uses 
the same notation as this paper. 

Let Bool = {0, 1}. As usual, let Z, M, C represent the set of integers (negative 
and non-negative), real numbers, and complex numbers, respectively. For integers a, 
b such that a < b, let Za^b = {ci, a + 1, . . .b—l,b}. For F equal to Z or M, let F^'^ and 
F-° represent the set of positive and non-negative F numbers, respectively. For any 
positive integer n and any set S, let denote the Cartesian product of n copies of 
S] i.e., the set of all n-tuples of elements of S. 

For any (not necessarily distinct) objects ai, 02, a^, . . ., let {ai, 02, a^, . . .}ord de- 
note an ordered set. For some object b, let b{ai, 02, a^, . . .}ord = {bai, 602, ba^, . . ■}ord- 
Let be the empty set. For an ordered set S", let S"^ be in reverse order. 

We will use 9{S) to represent the "truth function"; 9{S) equals 1 if statement 
S is true and if 5 is false. For example, the Kronecker delta function is defined by 
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5| = S{x, y) = 9{x = y). For x G Bool, 

1 

5^(-l)'^- = 25(x,0). (1) 

For any positive integer A^, we will use where z = l,2,...,A^to denote the 
standard basis vectors in dimensions; i.e., [cjjj = for i,j G Zi^tv- 

/„ and 0„ will represent the n-dimensional unit and zero matrices. 
For any matrix A and positive integer n, let 

A®" = A ® ■ • ■ ® A ® A , (2) 

n copies of A 

A®"" = A(B---(BA(BA . (3) 

71 copies of A 

For any matrix A G C™^" and p = 1,2, oo, || A||p will represent the p-norm of 
A, and its Frobenius norm. See ^I] for a discussion of matrix norms. 

Let X G C"^^. As is customary in the Physics literature, ||a;||2 will also be 
denoted by \x\ and called the magnitude of x. For any complex matrix A, we will use 
abs{A) to denote the matrix that is obtained from A by replacing each of its entries 
by its absolute value. In other words, [afes(y4)]jj = abs{Aij) = \Aij\. (Careful: Ref . [TT] 
and many other mathematical books call \A\ what we call abs{A)). 

Any a: G M can be expressed as a doubly infinite power series in powers of 
a base E G Z-^^: x = ± Yl'^=~oo ^aE". This expansion can be represented by: 
±(- ■ ■ a;iXo.2;_iX_2 ■ ■ ■ )\>e, which is called the base E representation of x. The plus or 
minus in these expressions is chosen to agree with the sign of x. It is customary to 
omit the subscript \)E when E = 10. For example 2.25 = 2 + | = (1.01)i,2 

Suppose X = (■ ■ ■ XiXo.a;_iX_2 ■ ■ ■ )b2 ^ R-^. Note that division by 2 shifts the 
binary representation of x one space to the right: | = {■ ■ ■ xi.xq ■ ■ ■)\,2- Likewise, 
multiplication by 2 shifts the binary representation of x one space to the left: 2x = 
(• • ■X-i.X-2 ■ ■ ■ )b2- In general, for any a e Z, ^ = {■ ■ -Xa-Xa-i ■ ■ ■ )b2- 

Define the action of an overline placed over an a G Bool by = 1 1 = 0. 
Call this bit negation. Define the action of an oplus placed between a,b E Bool 
by a © 6 = 6{a ^ b). Call this bit addition. One can extend the bit negation and 
bit addition operations so that they can act on non-negative reals. Suppose x = 
(■ ■ ■ xiXo.a;_ix_2 ■ ■ ■ )\,2, and y = {■ ■ ■ yiyo-y-iy~2 ■ ■ ■)\>2 are non-negative real numbers. 
Then define the action of an overline over x so that it acts on each bit individually; 
i.e., so that [ This overline operation is sometimes called bitwise negation. 

Likewise, define the action of an oplus placed between x and y hj (xQ) y)a = Xa® ya- 
This oplus operation is sometimes called bitwise addition (without carry). 

For any a; G M, the floor function is defined by \_x\ = max{j G Z : j < x}, 
and the ceiling function by [x] = min{j G Z : j > x}. For example, if x = 
(• ■ • xiXo.x_ix_2 ■ ■ ■ )b2, then [xj = (■ ■ ■ X2XiXo)b2- 
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We will often use Nb to denote a number of bits, and Ns = 2^^ to denote 
the corresponding number of states. We will use the sets Bool'^^ and Zo^at^-i inter- 
changeably, since any x G Z,o,Arg_i can be identified with its binary representation 
{xNg_i - ■ ■ XiXq\2 e Bool^^. 

For any x = {x^b-i ■ ■ ■ XiXo)\,2 ^ '^o,Ns-i, define x^ = {xqXi ■ ■ ■XNg_i)\,2; i.e., 
x^ is the result of reversing the binary representation of x. 

Suppose IT : Zo,Arg_i — > Zo,Ar<j_i is a 1-1 onto map. (We use the letter n to 
remind us it is a permutation; i.e., a 1-1 onto map from a finite set onto itself). One 
can define a permutation matrix M with entries given by My^ = 0{y = 7r(x)) for 
all x,y & Zo^AT^-i. (Recall that all permutation matrices M arise from permuting 
the columns of the unit matrix, and they satisfy M'^M = 1.) In this paper, we will 
often represent the map n and its corresponding matrix M by the same symbol vr. 
Whether the function or the matrix is being alluded to will be clear from the context. 
For example, suppose A is an Ns dimension matrix, and tt is a permutation on the 
set Z^ATg-i. Then, it is easy to check that for all i,j G Zo,Afs-i, {^'^A.)ij = and 

Suppose ttb '■ Zo^ATg-i — > Zq^^b-i is a 1-1 onto map (i.e., a bit permu- 
tation). TiB can be extended to a map ttb : Zq^^^^i — > Z^ATg-i as follows. If 
X = {xnb-1- ■ ■XiXo)\,2 e Zo,Arg_i, then let [7rB{x)]a = for all a G Zq^n^^i. 

The function ttb '■ '^o,Ns-i ~^ '^o,Ns-i is 1-1 onto, so it can be used to define a per- 
mutation matrix of the same name. Thus, the symbol tvb will be used to refer to 3 
different objects: a permutation on the set Zq^n^^i, a permutation on the set Zo,Ar<;_i, 
and an Ns dimensional permutation matrix. All permutations on Zq^Nb-i generate 
a permutation on Zo,Arg_i, but not all permutations on Zo,Arg_i have an underlying 
permutation on Zq^n^^i. 

An example of a bit permutation that will arise later is ttb', it maps nnii) = 
for all i G Zq^n^^i and HR^a) = Nb — 1 — a for all a G Zq^n^^i. 

3 Gray Code 

In this section, we will review some well known facts about Gray codejlS, . (Gray 
code was named after a person named Gray, not after the color.) 

For any positive integer Nb, we define a Grayish code to be a list of the 
elements of Bool^'^ such that adjacent A'^g-tuples of the list differ in only one com- 
ponent. In other words, a Grayish code is a 1-1 onto map Ticish '■ Z^at^-i — Zo,Arg_i 
such that, for all k G Zo^Ns-2, the binary representations of 7iGish{k) and 7iGish{k + 1) 
differ in only one component. For any Nb > 1, there are many functions TTash that 
satisfy this definition. 

Next we will define a particular Grayish code that we shall refer to as "the" 
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Gray code and denote by -kg- The Gray code for A^^ = 1, 2, 3 is: 









k 


{k),2 




k 


(A:)b2 


ncik) 







00 


00 













1 


01 


01 


1 


1 


1 




2 


10 


11 








3 


11 


10 



k 




1TG[k) 





000 


000 


1 
i 


UUi 


UUi 


2 


010 


Oil 


3 


oil 


010 


4 


100 


110 


5 


101 


111 


6 


110 


101 


7 


111 


100 



(4) 



The Gray code can be defined recursively as follows. Let Tq = 0. For A^^ > 0, let 
Tj\fg equal the set Bool^'^ ordered in the Gray code order. In other words, Vi 
{7rG(0),7rG(l),7rG(2), 



,7rG(2^«-l)W. Then 

= {or 



ord 



(5) 



for A''^ G Zo,oo-(See Section El for ordered set notation.) 

From the recursive definition of the Gray code, it is possible to prove that if 
k = {■ ■ ■ /c2fci^o)b2 and g = {■ ■ ■ g29igo)b2 are nonnegative integers such that g = ixaik)^ 
then 



ga — ka ® ka+l , (6a) 

for all a G Zq^oo- (For all a > Nb — 1, ka = ga = Eq.(|Ua|) specifies A''^ linear 
equations for the Nb components of g expressed in terms of the Nb components of 
k. These equations can be easily inverted using Gauss Elimination to get: 

ka = ga® ga+1 © © © " " " • (6b) 

Eqs.dni) can also be written in terms of the fioor function: 

k 

g = k(B[-\ , (7a) 

^ = ^7©Lfj©L|j©L|j©--- • (7b) 

As in Section 121 suppose tvb represents a permutation on Zq^at^-i which gener- 
ates a permutation on Zq^^s-i of the same name. Clearly, ttb o ttg is a Grayish code. 
Indeed, tib o ttg is a 1-1 onto map, and permuting bits the same way for all elements 
of a list preserves the property that adjacent A^s-tuples differ in only one component. 
(Note, however, that it is easy to find tt^'s such that ttg ° t^b is not a Grayish code. 
Hence, to preserve Grayishness, one must apply the bit permutation after 77^, not 
before). 
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4 Hadamard, Paley and Walsh Matrices 



In this section, we will review some well known facts about the so called Hadamard, 
Paley and Walsh matrices (a.k.a. transforms) JL2., . 

For any positive integer Nb-, we define the A'^g-bit Hadamard matrix by 



the iVe-bit Paley matrix by 



1 sr^Ng~l„Ng~l cJVs-1 



and the A'^B-bit Walsh matrix by 



where k,r E Z(^2^g_i, k = {■ ■ ■k2kiko)\,2 and r = (■ ■ ■r2riro)b2- We will often omit 
the subscript A''^ from H]\rg , P^r^ , W^g in contexts where doing this does not lead to 
confusion. 

Note that H, P, W are real symmetric matrices. 

For j G Zo,Arg-i, define the "reversal" function 'n'R{NB)U) = J^y the "nega- 
tion" function 71n{Nb)U) = J- The function 7VG(^j^g) for A^'^-bit Gray code has been 
defined previously. The functions ttr^Nb): '^n{Nb) T^GiNB) ^^^^ 1-1 onto so they can 
be used to define permutation matrices of the same name (See Section |21) We will 
often write tir, hn and tvg instead of vrij(Ar^), TT^^Ng) and tig{Nb) contexts where 
this does not lead to confusion. 

Note that tir and ti^ are symmetric matrices but ttc isn't. 

Next we will show that the A^^-bit Hadamard, Paley and Walsh matrices all 
have the same columns, except in different orders. More specifically, P, W are 
related to each other by the following equations: 

H'Kr = P, P'Kg = W. (11) 
A more pictorial way of expressing Eqs. ljllj) is: 

^i±Rpi±a^^ (12) 

Taking the transpose of both sides of Eqs. (fTT|) leads to 

nRH = P, TilP = W. (13) 

In the last equation, we have used the fact that matrices H, P, W, ttr, are sym- 
metric but ttg isn't. 
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and 



Comparing Eqs.(fTT|) and (fT^. we see that 

Htir = -krH , 



(14) 



In fact, Eq.lfT^ can be generalized as follows. Suppose tt^ is a bit permutation on 
^o,JVs-i- Then 



H.. 



1 



Hkr 



SO 



Httb = ttbH . 

Eq. (fT7jl becomes Eq. (fT^ when ttb = t^b.- 
To prove Eqs. fjllll . note that 
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ik = ^ Hij{'^R)jk 
3=0 



1 = 



Similarly, 



Pik ■ 



Ns-l 

PijMjk 

1 w » . . ^NjD—l 



1 



(16a) 
(16b) 

(16c) 
(16d) 

(17) 



:i8a) 



(18b) 



:i8c) 



(19a) 
(19b) 

(19c) 

(19d) 
(19e) 
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The square of H, P and W is one. Indeed, using Eq.dH) we get 



Ns-l 
j=0 



1 1 



2^3 ■•'X^X^'^ 1)^0.0 



(20a) 

(20b) 
(20c) 



and 



= i7TRH){H7TR) = 1 



= {7i^P){PnG) = 1 



(21) 



(22) 



Since their square equals one, and they are real symmetric matrices, H, P and W are 
also orthogonal matrices. 

From the definitions given above for ttrjTTnjT^g, one can prove by induction 
on Nb that these matrices obey the following recursive equations: 



T^R{Nb) ® (0, 1) 



1 

1 



and 



T^GiG) = 1 , T^G{Nb+1) 



T^GiNB) 

'T^G{Nb)'^N(Nb) 



(23a) 
(23b) 

(23c) 



Similarly, from the definitions given above for i/, P, VF, one can prove by induction 
on Nb that these matrices obey the following recursive equations: 



Np. 



1 1 



Fn = 1 , P, 



Pnb^{^A) 
Pn,^{1,-1] 



1 J 71' 
1 

71' 



and 



Wo = l, Wnb+i 



Wnb ® (1, 1) 



1 

71 



(24a) 
(24b) 

(24c) 
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By virtue of Eqs. (fTT|) . 



W = Httrttg , 



(25) 



Eq. (j25|) means that the permutation tt^ttg will permute the columns of H to give W. 
Expressing Eq. ()25|l in component form, we find 



Wij = Hir{-KR)rk{'^G)kj 
r,k 



r,k 



i,-KROTTc{j) ■ 



(26a) 

(26b) 
(26c) 



Thus, if we denote the columns of Hn^^ and Wnb by hj and Wj, respectively, then 



(27) 



for j G Zo,Ns-i- 



5 Constancy 

In this section, we will define a property of vectors called constancy. The columns of 
H^g can be conveniently classified according to their constancy. 
Consider the 3-bit Hadamard matrix: 
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(28) 



where we have labelled the columns of by /ij, where the index j is given in its 
binary representation. According to Eq.((27j), to get from Eq. ipHj) . one can simply 
reorder the columns of in bit-reversed Gray code. The columns of (and of W^) 
can be classified according to their constancy. We define the constancy C{h) of a 
vector h to be the smallest number of identical adjacent entries of h. For example. 
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C([l, —1, 1, 1]'^) = 1 and C([l, 1, —1, —1]'^) = 2. The next table gives the constancy of 
the columns of ifa, with the columns listed in the order in which they appear in W^. 
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TTr O TTnik) 
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001 
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(29) 



It is clear from Eq.(j2ni) that the columns of are listed in order of non-increasing 
constancy, and that the constancies of the columns of are all powers of 2. The 
literature on Walsh matrices often refers to the index that labels the columns of W 
as the sequency of that column. Thus, as sequency increases, constancy decreases 
or stays the same. Sequency and Constancy are analogous to Frequency and Period, 
respectively, in Fourier Analysis. 

Note that given any matrix A, more than one of the columns of A may have 
the same constancy. We will refer to: the number of columns of A with the same 
constancy K, as: the multiplicity of the constancy K in the matrix A, and 
denote it by M.a{C = K). In this paper, we are only concerned with the case where 
A equals the A^^-bit Hadamard matrix so we will henceforth omit the subscript A 
from M.a{C = K). Sometimes we will also abbreviate Ai{C = K) by A4{K), if doing 
this does not lead to confusion. The next table gives the multiplicity of the constancy 
K in the A^-bit Hadamard matrix: 



M{C = K) 


K = 1 


K = 2 


K = 4: 


K = 8 ■■■ 


Nb = 1 


1 


1 








Nb = 2 


2 


1 


1 





Nb = 3 


4 


2 


1 


1 



It is also convenient to define 

M{C>Ns')= Yl MK). (31) 

K>Ns' 

We shall call this the cumulative multiplicity of the constancy. The next table 
can be easily obtained from Eq. (j30|) and Eq. (j3T| . It gives A4{C > Ns') for the A^-bit 
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Hadamard matrix. 



M{C> Ns') 



Ns' 



2° Ns' = 2^ Ns' = 2'^ Ns' = 2^ 



Nb = 1 
Nb = 2 
Nb = 3 



2 
4 
8 



10 
2 10 
4 2 1 



(32) 



It is clear from Eq.^ that M{C > Ns') = 2^b-Nb'q(^^^ > ^^ere A^^' = 2^^'. 

6 Symmetries of Multiplexors 

In this section, we discuss some symmetries of exact decompositions of U (2)-multipIexors. 

For simplicity, we will first consider the i?jy(2)-multiplexors used in Ref. [5j. 
Ref. [Hj uses t/(2)-multiplexors that are more general than the i?y(2)-multiplexors 
used in Ref-ISj- At the end of the paper, we will discuss how to generalize our results 
for i?y(2)-multiplexors so that they apply to the more general multiplexors used in 



Below, we will present some quantum circuit diagrams. Besides the circuit 
notational conventions defined in Refs.jH] and (TU], the circuits below will use the 
following additional notation. A square gate with an angle 6 below the square will 
represent exp(z^(Ty) applied at that "wire". Typically, we will consider a SEO con- 
sisting of alternating one-qubit rotations and CNOTs. The SEO will always have a 
one-qubit rotation at one end and a CNOT at the other. The angle for the one-qubit 
rotation that either begins or ends the SEO will be denoted by ^oo...o- Given two 
adjacent angles 6^ and in the SEO, (6)b2 and {h')\,2 will differ only in one compo- 
nent, component a, where a is the position of the control bit of the CNOT that lies 
between the and 9y gates. 

If we take the Hermitian conjugate of the multiplexor exp(i Ylib&Booi^B-'^ 4'h(^y® 
Pb), and then we replace the angles (pb by their negatives (and also the angles 9b, 
Hadamard transforms of the (pb, by their negatives), we get the same multiplexor 
back. Henceforth, we will refer to this symmetry transformation as time reversal. 
Thus, an i?j^(2)-multiplexor is invariant under time reversal. 

Suppose hb is a bit permutation on A^^ — 1 bits. If we replace (pb by (pTrgib) 
(and also 9b by ^,rs(6)) and Pb by P^Bib) in the multiplexor exp(? ^^g^^^jjv^-i (pbay ® 
Pb), we get the same multiplexor back. Henceforth, we will refer to this symmetry 
transformation as bit permutation. Thus, an i?j^(2)-multiplexor is also invariant 
under bit permutation. 

FigO shows how time reversal and bit permutation act on a sequence of one- 
qubit rotations and CNOTs. More examples of the application of these transforma- 
tions will be given below. 



Ref.p. 
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TIME 
REVERSAL 



e e' 



e' 



□ — ^ 



bit a 



-a 



71, 



D- 



f bit TigCa) 







Figure 2: Examples of tlie action of time reversal and bit permutation on a string of 
one-qubit rotations and CNOTs. 



Figure 3: A half-moon node represents a projector Ph where b G Bool. The half- 
moon node may be omitted when it appears in a multiplexor whose t/( 2^^) -subset 
does not depend on the index b. This figure is an example of this principle. 



Recall from Ref.0 our definition of a general multiplexor with control 
qubits K and Nt target qubits r: J2b&Booi^K Ub{f)Ph{K). In a multiplexor whose 
matrices Uh are independent of the a component ba of 6, we can sum P^^ over b^ G 
Bool to get 1. Such a multiplexor acts as the identity on qubit a. When representing 
such a multiplexor in a circuit diagram, we can omit its half-moon node on qubit line 
a. FigOl shows a very special case of this principle, a special case that will be used 
in the circuit diagrams below. 

FiglU shows two possible ways of decomposing an i?j^(2)-multiplexor with one 
control. The decomposition (a) in FigEJis equivalent to: 



b^Bool 



(33) 



Let LHS and RHS stand for the left and right hand sides of Eq.()33|). Recall that 
n = diag{0, 1) = Pi and n = 1 — n = diag{l,0) = Pq. Eq. (jHH|) can be proven as 
follows: 
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01 Go 



(a) 



(b) 



Figure 4: Two possible decompositions of an i?j^(2)-multiplexor with 1 control. 



RHS 



j<Ty{i){eo+e^[Po{o)-Pim} 



= LHS . 

To arrive at Eq. (j34cp . we expressed 6*0, Oi in terms of (f>o, 0i using 
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' 1 1 " 




00 
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1 -1 





















(34a) 
(34b) 
(34c) 



(35) 



(1) 



(36) 



replace the angles 06 and 6b by their negatives, we get 

exp(z J2 My®Pb] =a,(l)"(°)e''^'^-«a,(l)"(°)e^^°'^- 

\ beBool J 

Eg ■ (j36j) is equivalent to decomposition (b) in FigEJ Thus, decompositions (a) and 
(b) in FiglU transform into each other under time reversal. 

FiglHl shows four possible ways of decomposing an i?y(2)-multiplexor with two 
controls. FiglHl was obtained by applying the results of FigsOl and El 

In Figini note that decompositions (a) and (b) transform into each other under 
time reversal. Decompositions (c) and (d) do too. Furthermore, decompositions (b) 
and (c) transform into each other under bit permutation. 

The decompositions exhibited in Fig^can also be expressed analytically. For 
example, decomposition (b) is equivalent to: 



exp i ^ 0fcO-y (g) Pb I 

\ beBool'2 J 



nil)A9io^yi2) 



n{Q)jeiiay(2) 



0^.(2) 



n{l)jeoiayi2) 



0^.(2) 



n(0) 



(37) 
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Figure 5: Four possible decompositions of an i?y(2)-multiplexor with 2 controls. 



Eq.(j27j) can be proven using the same techniques that were employed in Eqs. ljH^ to 
prove Eq.(j2Sl)- The proof requires that we assume: 
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(38) 



? — T T T T 

a 1 1 

a 1 

Figure 6: One of several possible decompositions of an i?j^(2)-multiplexor with 3 
controls. 
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Figure 7: One of several possible decompositions of an i?j^(2)-multiplexor with 4 
controls. 



Fig. ini (ditto, [7j) shows one of several possible decompositions of an Ry{2)- 
multiplexor with 3 (ditto, 4) controls. In general, decompositions for multiplexors 
with Nk controls can be obtained starting from decompositions for multiplexors with 
Nk — 1 controls. 

7 Approximation of Multiplexors 

In this section, we finally define our approximation of multiplexors. We give the 
number of CNOTs required to express the approximant, and an upper bound to the 
error incurred by using it. 

So far we have used A''^ to denote a number of bits, and Ns = 2^^ to denote 
the corresponding number of states. Below, we will use two other numbers of bits, tje 
and ris', where tje = Nb — 1 and t]b' < rjB- Their corresponding numbers of states 
will be denoted by r]s = 2''^ and rjs' = 2''^ . 

Define an r]s dimensional matrix V by 

V = Httb-kg , (39) 

where ttb is an arbitrary bit permutation on rjB bits. Eq. (j39|) is a generalization of 
Eq.(j2Hl). Both equations define a new matrix (either V or the Walsh matrix W) by 
permuting the columns of the Hadamard matrix H. V becomes W if we specialize 
the bit permutation tib to tir. If we denote the columns of V by Vj for j G Zo^rjs-iy 
then 



which is the counterpart of Eq. ()27|) . 

In Ref. [5 , the decomposition 
following Hadamard transform: 

e = 



— hirgOTTaij) 5 (40) 

of an Ry{2) multiplexor starts by taking the 
H,J. (41) 
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where r]B = Nb — 1 and -qs = 2^'^ . The vectors {wjjvj constitute an orthonormal basis 
for the space M^-^ in which (p lives, so (p can always be expanded in terms of them: 



'7S-1 
i=0 



(42) 



Now suppose that we truncate this expansion, keeping only the first rjs' terms, where 
rjg' = 2^^ and -qs' G Zq^tva-i- Let us call 0' the resulting approximation to (j): 



Define 9', an approximation to 9, as follows: 

1 



(43) 



9' 



(44) 



If we let {ej}vi denote the standard basis vectors, then 

hi 

h\ h 



H V- 



■^BOT^cii) ^■n-B°'^G{i) 



(45) 



Therefore, 



9' 



°TrBOTTc{i) \ 



(46) 



By virtue of Eq.pUj). if we list the components {^[,}v6 of 9' in the Grayish 
code order specified by the map ttb o ttq, then the items in the list at positions 
from r]s' to the end of the list are zero. Consider, for example, FiglHl which gives 
the exact decompositions for a multiplexor with 2 controls. Suppose that in one of 
those decompositions, the angles 9b s in the second half (i.e., the half that does not 
contain 6'oo) of the decomposition are all zero. Then the one-qubit rotations in the 
second half of the decomposition become the identity. Then the three CNOTs in the 
second half of the decomposition cancel each other in pairs except for one CNOT that 
survives. The net effect is that the decomposition for a multiplexor with 2 controls 
degenerates into a decomposition for a multiplexor with only 1 control. The number 
of control bits is reduced by one in this example. In general, we can approximate 
a f/(2)-multiplexor by another t/(2)-multiplexor (the "approximant" ) that has fewer 
controls, and, therefore, is expressible with fewer CNOTs. We will call the reduction 
in the number of control bits the bit deficit 6b- Hence, Sb = r]B — Vb'- 
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If NcNOT denotes the number of CNOTs in an approximant with bit deficit 
6b, then it is clear from FigsjH ElandlTlthat: 



Sb = 





1 


2 




Nb-2 


Nb-1 


NcNOT = 


2^3-1 


2^3-2 


2^3-3 




2 






(47) 



Hence, for 5b G 'Lo,Nb-2, Ncnot = 2^s-i-<5s^ but for 6b = Nb - 1, Ncnot = 0. 

The bit permutation tt^ on which the approximation of a multiplexor depends 
can be chosen according to various criteria. If we choose hb = ^r, then our approxi- 
mation will keep only the higher constancy components of <p. Such a smoothing, high 
constancies approximation might be useful for some tasks. Similarly, if we choose 
vtb = 1, then our approximation will keep only the lower constancy components of 
0, giving a low constancies approximation. Alternatively, we could use for ttb a 
bit permutation, out of all possible bit permutations on r/^ bits, that minimizes the 
distance between the original multiplexor and its approximant. Such a dominant 
constancies approximation is useful if our goal is to minimize the error incurred 
by the approximation. 

The error incurred by approximating a multiplexor can be bounded above as 
follows. Let {e^'^'''^^}\/b£Booi'iB denote the i?j^(2)-subset of an i?j^(2)-multiplexor T and 
{e'''^b'^y}\/heBooi''B that of its approximant T'. Call ||T' — T||2 the error of approximating 
T by T'. Note that 



||T' - T||2 = II Qb^Booi^B {e^^'y - e''*>^''y)h (48a) 

= maxlle^'^^'"^ -e^'^^^HU (48b) 
b 

< max \<f)'b - (f)b\ = U' - 0||oo . (48c) 

b 

To arrive at step Eq. ()48c|) . we used the results of Appendix |X1 We will sometimes 
refer to ||0' — 0||oo as the linearized error, to distinguish it from the error ||T' — T||2. 

A simple picture emerges from all this. The error e and the number of CNOTs 
Ncnot are two costs that we would like to minimize. These two costs are fungible 
to a certain extent. Given a multiplexor T, and an upper bound eo on e, we can use 
Eqs.flTTj) and (PH|) to find the approximant T' with the smallest Ncnot- Similarly, 
given a multiplexor T, and an upper bound {Ncnot)o on Ncnot, we can use Eqs. ljlTj) 
and (j48p to find the approximant T' with the smallest e. 

At this point we encourage the reader to read Appendix ^ It discusses the 
output of a computer program that calculates 0' from (p via Eq. ()43|) . 

Next we will show that Eq. pHjl can be simplified considerably by taking into 
account the explicit values of the column vectors vj. 

To get a quick ghmpse of the simplification we seek, consider first the special 
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case rjB = 2. We have 



H2 = - 
2 



^00 


hoi 




hii 








hii 


hoi 


+1 


+1 


+ 1 


+ 1 




+ 1 


+ 1 


+ 1 


+ 1 


+1 


-1 


+ 1 


-1 


+ 1 


+ 1 


-1 


-1 


+1 


+1 


-1 


-1 




+ 1 


-1 


-1 


+ 1 


+1 


-1 


-1 


+ 1 




+ 1 


-1 


+ 1 


-1 



(49) 



Define a matrix jj, by 



1 1 
1 1 



(50) 



For any matrix A, let A{:, i : j) be the submatrix of A obtained by keeping only its 
columns from i to j. It is easy to check that 



W2(:,0 : 3)W2(:,0 : 3)' = diag{l, 1,1,1 



and 



W2{:,0:1)W2{:,0:1) 



W2{:,0: 0)PF2(:,0 : 0)'^ 



T 



/i 




-(/i © /i) 



(51) 
(52) 

(53) 



In each case, we formed a "decimated matrix" W{:, : rjs' — 1) from where rjs' = 
2VB _ Then we showed that the projection operator W{:, : rjs' — 1)W{:, : r]s' — 1)^ 
onto the column space of the decimated matrix, is a matrix whose entries are all 
either or ^jj, and these entries sum to one along each row (or column). Given a set 
S of real numbers, and given 5*1 C S, call the average of the elements of 5*1 a "partial 
average" of the elements of S. For example, if rjs' = 2 and ttb = t^r, then V = W and 
(j)' = W2{:, : 1)W2{-, : l)'^(f). From Eq.(jS2), the entries of (p' are partial averages of 
the entries of (p. 

Next we show how to simplify Eg . (|43|) for arbitrary rjB, not just for rjB = 2. 
Define an 7^5 dimensional matrix F^''^'^^'-' by 



r=0 



(54) 



Below, we will show that T^^^'^^''^ reduces to: 



qs 



25b J L2<5sJ) ■ 



(55) 
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For example, when ijb = '2,r]B' = 1, Eq.(|^ becomes 



[r?/^]=^[^(L|j = L^J) 
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00 
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01 
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11 








1 


1 



(56) 



Eqs. fl54|) and are given in component form. The identical statements 
written in matrix form are, respectively. 



W, 



VB 







vs-vs' 



T 

Vb ' 



and 



(57a) 
(57b) 



Y(Vb,Vb') = ^ - 
2<5s 



\3VS 



Eq. ()55p can be proven from Eq. ()54|) as follows: 



(5^ 



-J- 

Y(vb,Vb') = ± ^ (_l)E:£o'E^2o'('?.r^+r-,«.){5::5-V<5:?^) (59^) 
r=0 

if' ' '\ 

= —\ Yl ■ ■ ■ 5Z 5Z ^(''VB-I = = ■ ■ ■ = r^B-5B = 0) 



r,,^_i=0 ri=Oro=0, 

(_l)Z.«=o 2^,3=0 '•/3ta(,6^+^ (where t = q ® s) (59b) 

= ~\ ' ^(''vB-i = = ■■■ = rri^-Ss = 0) 

\r,jg-i=0 ■n=Oro=oJ 

(_l)^')B-l(*0+tl)+»-r;s-2(tl+t2) + ---ri(t^^_2+t,,g-l)+ro(t,,g-l) (59c) 

— ]_onB-SB s:'''^B-^ sj^'vB-^ . . . C^Qd) 

= ^0(1^1 = 1^1) ■ (59e) 

To arrive at Eq. ()59d|) . we used Eq.((T)). 

Recall that W and V can both we obtained by permuting the columns of H: 

W = Htirttg , (60) 
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and 



V = HnsTTG ■ (61) 

From these two equations and from the fact, proven earher, that tib and tir commute 
with H, we get 

V = -kb-krH-kr-kg = ttbtj^rW . (62) 

Thus, 



VB 







= -TTB-^RWrig 

Hence, if we define if and tf' by 



ns-vs 

Ins' 
0, 



vs-vs 



W^„'^r7tI$. 



(63a) 
(63b) 



then 



^ = TTRTT^cj) , tp' = 7lR7rB<P' 



(64) 



(65) 



We see from Eq. (l65|l that even when tir 7^ t^r, the entries of i/j' (which are the 
same as the entries of (p' but in a different order) are partial averages of the entries 
of ip (which are the same as those of but in a different order). Thus 



min( 

k 



h) < 



A' < max(0fc) 

k 



for all j G TiQ^ris-i- This last equation implies 



10',- -0i 



< max( 

k 



mm 

k 



(66) 



(67) 



for all j. 

As we mentioned before, the quantum compiling algorithm of Ref. [6J uses U (2)- 
multiplexors that are more general than the i?y(2)-multiplexors considered above. 
Luckily, the above results for i2j,(2)-multiplexors are still valid, with minor modifi- 
cations, for the more general ones. Indeed, the [/(2)-subset of the multiplexors used 
in Ref.|ni is of the form {exp[i{(f>if,asi + <p2bO's2)](^w''^}yb- Ref.|B3 defines vectors 0i 
and 02 from the parameters {(f)ib}\/b and {02fe}vfe5 respectively. It then defines 9i 
and 02 as Hadamard transforms of (pi and (p2, respectively, just as Eq. ()41|) defines 
^ as a Hadamard transform of (p. We can define approximations (p'l and 9[ by re- 
placing (p, 9, (p', 9' by 01, 9i, (p[, 9[, respectively, within Eqs.(PH|) and (jUj). We 
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can define approximations 02 ^2 analogously. The expansions of (f)\ and (p'2 in 
the Vi basis can be truncated at the same ijs'. The table given in Eq.(|T7j) for the 
number of CNOTs still applies, except that Ncnot may change by 1 if we elim- 
inate the ^o...oo gate as in Ref.|ni. When a [/(2)-multiplexor T with ?7(2)-subset 
{exp[i(0ibcrsi + 4'2bO's2)Ww''^}vb is approximated by a [/(2)-multiplexor T' with U{2)- 
subset {exp[i{(p[i^asi + (f)2h(^s2)]'^w^^}vb, one can show, using the results of Appendix 
m that 



IT' — Til 2 < max 



(68) 



which is a generalization of Eq. 



A Appendix: Distance between 
two SU{2) matrices 

In this appendix, we establish a well known(see Ref.jTT], page 574) upper bound for 
the distance (measured in either the 2-norm or the Frobenius norm) between two 
SU{2) matrices. 

Let a, a' be 3d real vectors. Define Aa = a' — a. If | Aa| << 1, then 



= ||e^°'-V^°-^- III2 (69a) 
^ \\iAa-a\\2 (69b) 
= |Aa| . (69c) 

Next, we will show that this approximation can be turned into an inequality. 

Consider first the special case where a and a' both point in the Y direction. 

Then 



e^^""^-l||2 (70a) 

cos(Aa) — 1 sin(Aa;) 
— sin(Aa) cos(Aa) — 1 



(70b) 



. Aa / -sin(^) cos(^) , . 
|2sm(— ) \L\ .;'/aA II2 (70c) 



2 ' \ -cos(^) -sm( 2 

= 2|sin(^)| (70d) 
< |Aa| . (70e) 
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To find an upper bound for ||e*°'' '^ — e'" '^||2 wlien eitlier a or a' does not point 
in tfie Y direction, we will use tlie following identity. For A, E ^ (Qnxn ^ ^ M-'^, 

To prove Eq. lffTjl . let C and TZ stand for the left and right hand sides of Eq. (|7T|l . It 
is easy to verify that 

(/:-7^)(0) = 0, ^(^ - ^) ^ _ 7^) ^ ^72) 

This initial value problem has the unique solution C — TZ = 0. 

In Eq. fj7ip . set A = ia-a, and E = iAa ■ a, and take the 2-norm of both sides. 
This yields 



e ' — e 2 



< / rfs ||e^"-"(^-^)||2|UAa-(T||2||e 



i{d+Aa)-(7S I 



lAal . 



(73a) 
(73b) 



One can also find an upper bound for the distance, in the Frobenius norm, 
between two SU{2) matrices. If U = e*° '^ and U' = e*"' "", then the eigenvalues of 
U are e*^, e~*^, for some real number 6. Likewise, the eigenvalues of U' are e*^ , e~*^ . 
Thus ti^U' — U) is real. If we denote the eigenvalues oi U' — U hy x^iy with ?/ G M, 
then ([/' — U)'^{U' — U) has a single eigenvalue + y"^ with algebraic multiplicity 2. 
Thus 

_ ei"-^||^ = ^ II gi"'-'? _ e*"-^||2 . (74) 

But we've already proven that He*"' "" — e*"'°"||2 is bounded above by |Aa| so 



B Appendix: Computer Results 

In this appendix, we discuss a simple computer program that verifies and illustrates 
many of the results of this paper. Our program is written in the Octave language. 
Octave is a gratis, open-source interpreter that understands a subset of the Matlab 
language. Hence, our program should also run in a Matlab environment with few or 
no modifications. 

Our main m-file is called my_moo.ni. When you run my_moo. Octave produces 
two output files called out_phis.txt and out_error.txt. 
A typical out_phis.txt file reads: 
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phi(l)= 0.133765891 
phi (2)= 0.270447403 
phi (3)= 0.307625920 
phi (4)= 0.311291575 
phi (5)= 0.452735037 
phi (6)= 0.569045961 
phi (7)= 0.653136015 
phi (8)= 0.867156088 

permutation 1 = (1,2,3) 
delta_B, phi_prime= 

0.134 0.270 0.308 0.311 0.453 0.569 0.653 0.867 

1 0.293 0.420 0.480 0.589 0.293 0.420 0.480 0.589 

2 0.387 0.504 0.387 0.504 0.387 0.504 0.387 0.504 

3 0.446 0.446 0.446 0.446 0.446 0.446 0.446 0.446 



permutation 2 = (1,3,2) 
delta_B, phi_prime= 

0.134 0.270 0.308 0.311 0.453 0.569 0.653 0.867 

1 0.293 0.420 0.480 0.589 0.293 0.420 0.480 0.589 

2 0.356 0.356 0.535 0.535 0.356 0.356 0.535 0.535 

3 0.446 0.446 0.446 0.446 0.446 0.446 0.446 0.446 



permutation 3 = (2,1,3) 
delta_B, phi_prime= 
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permutation 4 = (2,3,1) 
delta_B, phi_prime= 

0.134 0.270 0.308 0.311 0.453 0.569 0.653 0.867 

1 0.202 0.202 0.309 0.309 0.511 0.511 0.760 0.760 

2 0.356 0.356 0.535 0.535 0.356 0.356 0.535 0.535 

3 0.446 0.446 0.446 0.446 0.446 0.446 0.446 0.446 



permutation 5 = (3,1,2) 
delta.B, phi_prime= 

0.134 0.270 0.308 0.311 0.453 0.569 0.653 0.867 

1 0.221 0.291 0.221 0.291 0.553 0.718 0.553 0.718 

2 0.256 0.256 0.256 0.256 0.636 0.636 0.636 0.636 

3 0.446 0.446 0.446 0.446 0.446 0.446 0.446 0.446 



permutation 6 = (3,2,1) 
delta_B, phi_prime= 

0.134 0.270 0.308 0.311 0.453 0.569 0.653 0.867 

1 0.202 0.202 0.309 0.309 0.511 0.511 0.760 0.760 

2 0.256 0.256 0.256 0.256 0.636 0.636 0.636 0.636 

3 0.446 0.446 0.446 0.446 0.446 0.446 0.446 0.446 

The corresponding out_error.txt file reads 

error as function of (permutation\delta_B) 

12 3 

1 l.llOe-16 2.779e-01 3.627e-01 4.215e-01 

2 l.llOe-16 2.779e-01 3.324e-01 4.215e-01 

3 l.llOe-16 1.491e-01 3.627e-01 4.215e-01 

4 l.llOe-16 1.070e-01 3.324e-01 4.215e-01 

5 5.551e-17 1.491e-01 2.316e-01 4.215e-01 

6 5.551e-17 1.070e-01 2.316e-01 4.215e-01 

In this example, A^^ = 4 so r/^ = 3 and rjs = 8. The first 8 lines of 
out_phis.txt give the components of (f). In this case, the computer picked 8 in- 
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dependent random numbers from the unit interval, and then it sorted them in non- 
decreasing order. my_moo.m can be easily modified so as to allow the user himself to 
supply the components of (p. 

After listing ^, out_phis.txt lists the r/^! permutations tib of r/s bits. For 
each TiBi it prints the components of 0', listed as a row, for each value of 5b(=i'ow 
label). Note that for 6b = 0, (p' = 0, and for 6b = rjB, all (p'j are equal to the 

average of all the components of 0. Note also that for all values of 6b and j, (p'j G 
[minfc(0fc),maxfc(0fc)]. 

The second output file, out_error.txt, gives a table of the linearized error 
\\4>' " 0l|oo as a function of permutation number(=row label) and (5B(=column label). 
As expected, the error is zero when 6b is zero, and it is independent of the permutation 
7Tb when 6b is maximum (When the bit deficit 6b is maximum, the approximant has 
no control bits, so permuting bits at positions Zo^rjg-i does not affect the error.) 

Note that in the above example, the last permutation minimizes the error for 
all 6b- This last permutation is ttb = ttr = (bit-reversal), and it gives a high constan- 
cies expansion. Recall that for this example, my_moo.m generated iid (independent, 
identically distributed) numbers for the components of (p, and then it rearranged them 
in monotonic order. When (f) is chosen in this way, the graph {{j, (f>j)}vj has a high 
probability of lying close to a straight line, and a high constancy staircase is the best 
fit for a straight line. For this reason, almost every time that my_moo .m is operated in 
the mode which generates iid numbers for the components of 0, the high constancies 
expansion minimizes the error for all 6b- However, this need not always occur, as the 
following counterexample shows. Try running my_moo.m for Nb = 5, and for with 
its first 7 components equal to and its 9 subsequent components equal to 1. For 
this 0, and for 6b = 3, the high constancies expansion yields an error of 7/8 while 
some of the other expansions yield errors as low as 5/8. 

Note that although my_moo.m visits all t/^! permutations of the control bits, 
visiting all permutations is a very inefficient way of finding the minimum error. In 
fact, the r]B\ control bit permutations can be grouped into equivalence classes, such 
that all permutations in a class give the same error. It's clear from FigQthat we only 
have to visit (^^) = ^Jj^^,, (recall r]B' = i]b — 6b) equivalence classes of permutations. 
Whereas r^^! ~ r]B^'^ = e^'^^^^^ is exponential in r]B-, (5^) is polynomial in rjB for two 
very important extremes. Namely, when 6b or rjB is of order one whereas rjB is very 
large. Indeed, if 5^ = 1 or 77^' = 1, then (^^) = rjB] if (5^ = 2 or 77^' = 2, then 

iZ) = '-^^ etc. 
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